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Abstract 

The fundamental ingredients of the MCAS (multi-channel algebraic scattering) 
method are discussed. The main feature, namely the application of the sturmian 
theory for nucleon-nucleus scattering, allows solution of the scattering problem 
given the phenomenological ingredients necessary for the description of weakly- 
bound (or particle-unstable) light nuclear systems. Currently, to describe these 
systems, we use a macroscopic, collective model. Analyses show that the couplings 
to low-energy collective-core excitations are fundamental but they are physically 
meaningful only if the constraints introduced by the Pauli principle are taken into 
account. For this we introduce in the nucleon-nucleus system the Orthogonalizing 
Pseudo-Potential formalism, extended to collective excitations of the core. The 
formalism leads one to discuss a new concept, Pauli hindrance, which appears to 
be important especially to understand the structure of weakly-bound and unbound 
systems. 



1 Sturmians 

We outline here the use of Sturmian functions in the formulation of the Coupled-Channel 
scattering problem. Sturmians provide the solution of the scattering problems by matrix 
manipulation (hence the word "algebraic" in MCAS). They are known and used also 
in atomic and molecular physics, chemistry and field theory [TJ [2] . They provide an 
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efficient formalism for determination of S-matrices, scattering wave functions, bound 
states, and resonances. They work well with non-local interactions (such as those 
non localities arising from the effects due to Pauli exchanges). They allow a consistent 
treatment of Coulomb plus nuclear interactions, as well as the inclusion in the scattering 
process of coupled-channel dynamics. This occurs, for instance, when strong coupling 
to low-lying excitations of the target nucleus have to be taken into account. 

Sturmians (also known as Weinberg states) represent an alternative way to formulate 
the Quantum Mechanical problem. 

Consider a two-body like Hamiltonian H = H Q + V: then the Schrodinger equation 
is written in the standard (time-independent) way 

[E- H )*e = V^ E , (1) 

where E is the spectral variable, and f ^ is the eigenstate. 
Sturmians, instead, are the cigcnsolutions of: 

(E - H )^(E) = -^-^(E) , (2) 
Vi\E) 

where E is a parameter. The eigenvalue r\i is the potential scale. Thus the spectrum 
consists of all the potential rescalings that give solution to that equation, for given 
energy E, and with well-defined boundary conditions. 

The standard boundary conditions for Sturmians $i(E) are: 

E < Bound-state like; normalizable. 

E > Purely outgoing/radiating waves; non-normalizable. 

The spectrum of eigenvalues is purely discrete, and bounded absolutely. For short- 
range (nuclear-type) potentials, the eigenvalues can accumulate around only. 
Then, the single-channel S*-matrix can be written as 

S{E) - w (3) 

Alternatively, introducing the factor Xi(E,k) in momentum space 

Xi(E,k)=<k,c\V\$ i (E)>, (4) 
the 5-matrix can be rewritten also as 

S(E) = 1 - ink £ Xl ; k) 1 _^ {E(+)) Xi ; *) (5) 



Most interestingly, the last expression can be generalized to coupled-channel dynam- 
ics: one starts from a coupled-channel Hamiltonian with potential V cc i , and obtains an 
5-matrix of the form [3] : 

S cc ,(E)=S cc >-m^/k^k F c y / Xci(E^-,k c )- L^^(£ (+) ;fc c ') (6) 

i l-ViiE^ 1 ) 

It is remarkable that the following interpretation can be given to the last expression: 
the scattering process initiated in the asymptotic channel c is "captured" into Sturmi- 
ans. Subsequently the Sturmian propagates freely in the interaction region and finally 
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decays into the outgoing channels. This structure naturally reflects the description of 
the scattering process in terms of compound nucleus formation, and leads to an expres- 
sion which is rather similar to that obtained in i?-matrix formalism. However, in the 
Sturmian approach, such structure emerges directly from the Hamiltonian, while in the 
R-matrix formalism, the resonant (compound) structure is modeled phenomenologically 
in terms of specific boundary conditions given at the surface of an hypothetic R-space 
sphere. 

1.1 Resonances and bound states in terms of Sturmians 

Resonant structures as well as bound states can be obtained in terms of the properties 
of Sturmian eigenvalues. A bound state occurs when one of the eigenvalues moves 
toward the right on the real axis, and crosses the value 1 at some negative energy. That 
particular energy value corresponds to the bound-state energy. A resonance occurs 
when the eigenvalue, initially progressing along the the real axis, becomes complex 
before reaching the point (1,0). Such occurs for positive energies as the scattering 
threshold (E = 0) must be passed. The energy centroid of the resonance is the energy 
corresponding to the real part of sturmian eigenvalue matching the value 1. The width 
of the resonance can also be determined geometrically by the patterns of the sturmian 
trajectories, and relates to the imaginary part of the sturmian eigenvalue at the resonant 
energy. Such patterns can be seen in Fig. Q] In Fig. [2] one observes the situation in 
a realistic case, namely n- 12 C elastic scattering. The eigenvalue trajectories produce 
two 3/2 + resonances that can be clearly seen in the experimental data as well as in the 
theoretical calculation. 




Figure 1: How resonances and bound states are found in Sturmian theory. 



2 Model Coupled-channel potential and OPP 

To date, we have used a macroscopic potential approach: a nuclcon is scattered by a 
nucleus (light nuclei with + g.s. are considered) and we include couplings to first core 
excitations of collective nature (quadrupole, octupole, etc), since these couplings play 
an important role in the dynamics. The coupled-channel potential that describes the 
dynamics (including the couplings to collective-type target excitations) is an expansion 
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Figure 2: Neutron- 12 C in the 3/2+ channel: a realistic case. Low-energy resonances in 
3/2+ n- 12 C system. Sturmian patterns (left) and 3/2 + resonant cross-section (right). 
The dashed line in the left panel denotes the unit circle. 



over four operators, where the first two play a dominant role (central and spin-orbit), 
while the remaining two operators (orbit-orbit and spin-spin) produce small phenomeno- 
logical corrections. The sum over the various operator forms can be written as follows [4] 

VcAr) = y « < (^ J ; jn \°nfn(r, R, e itR )\(l's)j'l'; J* > 

n=C,ls,U,sI 

For all operators, the functional forms are expanded to second order in the core- 
deformation parameter. (For simplicity, we discuss here the case of a single quadrupole 
deformation effect: R = Rq(1 + fel%(d).) 

f n (r, R, 9) = /(°)(r) - (3 2 R P 2 (9)^fi°\r) 
The radial forms /^°^ (r) are spherically-symmetric functions. 

For the central, orbit-orbit, and spin-spin terms we consider a standard Wood-Saxon 
form: 

/W(r) - [1 + exp^]- 1 for (n = C,ll,sI) (7) 

For the spin-orbit term we consider Oi s —Is and not the full Thomas term, with 
the following radial form: 

/i°2w = ix[ 1 + ex p £ ^r 1 f° r a«) (8) 



2.1 First application: n- 12 C 

In our first application, we considered the scattering of neutrons off 12 C coupling the 
ground state of the target to the first two low-lying excitations 2^~ (4.43 MeV) and 0^ 
(7.63 MeV), and searched parameters to obtain a description of the resonant spectra 
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and scattering cross sections [H [5J However, many deeply-bound spurious states 
occurred in the bound spectrum. It was not possible to obtain a consistent description of 
both bound structure and scattering data with the same Coupled-Channel Hamiltonian. 
These spurious deeply-bound states originate from the violation of the Pauli principle. 
The phase-space corresponding to target nucleons in fully occupied shells has to be 
inhibited to the incoming nucleon. In the original CC model such condition was missing. 

Various methods have been suggested to remove the deeply-bound forbidden states. 
A recent article [7] on the subject contains an historical review on the various ap- 
proaches and their connections. In phenomcnological macroscopic-type calculations, 
such Pauli condition has been implemented first by introducing the Orthogonality Con- 
dition Model [5]. Alternatively, the orthogonality condition can be introduced directly 
in the Hamiltonian by addition of a new term in the potential, the highly non-local 
Orthogonalizing Pseudo-Potential [9]. With the advent of super-symmetric quantum 
mechanics [10], it was possible to define super-symmetric transformations [11] that pro- 
duce new local (and highly singular) potentials which also generate spectra free of 
spurious states. 

In our MCAS approach, we use the technique of Orthogonalizing Pseudo Poten- 
tials [5], which eliminates the deep bound states adding a new term in the nuclear 
potential. 

The "complete" nuclear potential for the n + 12 C case we considered, has the form 
(in partial-wave decomposition) 

V cc ^r,rO = F cc /(r)5(r-rO + <5cc'A c A c ^ 

A c (r) are the Pauli- forbidden deep (CC-uncoupled) bound states. The A parameter 
eliminates the deeply-bound spurious states, since a state in the OPP approach is 
forbidden in the limit A — > +oo, while it is allowed when A — > 

Thus, in 13 C, we considered two shells (the Osi and 0p| shells) to be Pauli blocked. 

2.2 The A-dependence 

One important aspect in the OPP method it to assess the behavior with respect to 
the A parameters. Originally, A was set to 100 MeV since this was sufficient to remove 
all the spurious states, but later it was found that some resonances were still sensitive 
to higher values of A, corresponding to a stronger orthogonality condition, and the 
energy centroids of those selected resonances improved by using larger values of A. In 
particular, a narrow 5/2~ resonance, which lies very close to threshold, illustrates this 
well. When A = 100 MeV, it is almost a zero-energy bound state. For higher values 
of A that resonance converges to a position in agreement with the peak observed in 
the evaluated data around 0.1 MeV. Also the energy centroid of the 1/2 - state around 
2.9 MeV stabilizes with higher A values. These effects indicate of the need of strong 
orthogonality conditions in the system. For further details, including the parameters 
set for the potential, we refer the reader to Refs. [HE!]. 

3 Applications of MCAS theory 
3.1 Analyzing powers of nucleons off 12 C 

Analyses of n- 12 C system were published in Refs. [UG3[B]. Note that all the parameters, 
including the spin-orbit term, have been fixed on the known spectrum of 13 C. Shortly 
after our analyses have been completed, new spin-polarization data measured at TUNL 
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Figure 3: The A dependence (of the elastic cross section for n- 12 C scattering) with 
values 100, 300, 500, and 1000 Mev (top-left, top-right, bottom-left and bottom-right, 
respectively) 



V (±) = -45.0 MeV 


Vjf ] = 0.42 MeV 




= 7.0 MeV 


#0 = 3.1 fm 


ag = 0.65 fm 


p2 = 


-0.50 



Table 1: Coupled-channel parameters for n- 14 C/p- 14 



were published |12j . Without any parameter adjustments or tuning, we could reproduce 
those low-energy data, as shown in Fig. dJ Thus, this was a test of our spin-orbit 
parametrization and a validation of the spin-structure of our interaction [13] . 

3.2 Low-lying unbound states of 15 C and 15 F 

Our second analysis concerned an unbound nucleus, 15 _F, whose properties were studied 
in connection with its weakly-bound mirror partner 15 C. Our analysis was triggered by 
recent low-energy experimental data on the 14 0-proton system, which we could analyze 
by using inverse kinematics. We included in the model low-lying excitations of 14 0/ 14 C, 
in a macroscopic coupled-channel model with parameters given in Table [1] 

Our results are shown in Fig.O In the evaluations, we used known properties of 15 C 
to predict new states in 15 F, in particular three narrow resonances of negative parity 
(l/2~, 5/2~ and 3/2~ in the range 5-8 MeV). To obtain these results, a new concept 
had to be introduced [2]: it is the concept of Pauli- hindered states. Up to now we 
have considered states that are Pauli allowed (A ~ 0) or Pauli forbidden (A > lGeV). 
Now we introduce a state that is neither prohibited nor allowed but simply suppressed, 
or hindered, by the Pauli principle. We expect that this situation can apply in weakly- 



6 




Figure 4: Neutron analyzing power A y at various lab energies in the MeV ranges (a) 
1.9, (b) 2.2, (c) 2.8, (d) 3.2 for the left panel, and (a) 3.41, (b) 3.62, (c) 3.78, (d) 3.92 
for the right panel. 




1 2 3 4 5 
E c m (MeV) 



Figure 5: Left: 15 F (top, p - 14 O differential cross-section 9 cm = 180°) & 15 C (bottom, 
the bound spectrum and the elastic scattering cross-section). Right: The resonant 
spectra for 15 C and 15 F 
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bound/unbound light systems where the formation of shells may become critical. These 
hindered states can be conveniently described in the OPP scheme with values of A of a 
few MeV. Thus the scheme of phase-space that is accessible to the proton interacting 
with the 14 core is: 

Pauli Forbidden (A ~ IGeV): 

0si /2 +0+ 0s 1/2 +0+ 0s 1/2 + 2+ 

0P3/2 + Of 0p 3 /2 + 0+ lf> 3 /2 + 2t 

qpi/2 + oi 

Pauli Hindered (A ~ 1 - 50MeV): 



Qpi/a + OJ Qpi/a + 2^ 

Pauli Allowed (A = OMeV): 

lsi/2+0^ ls 1/2 +0t lsi /2 +2+ 

04/2 + 0t 0d 5/2 + 0+ 0d 5/2 + 2+ 

0d 3/2 +0^ 0d 3/2 +0+ 0d 3/2 +2+ 



The need to consider intermediate situations between Pauli blocking and Pauli al- 
lowance has been registered before in the literature, mostly in connection with RGM 
approaches [151 US] . 

3.3 P-shells in mass=7 nuclei 

As a third application, we consider the spectral structure the mass=7 isobars: 7 He, 7 Li, 
7 Be, and 7 B. We describe these systems in terms of a single nucleon-nucleus interaction, 
explicitly including the low-lying core excitations of the mass-6 sub-system. Thus, we 
use the MCAS approach to determine the bound and resonant (above nucleon emission) 
spectra starting from a collective coupled-channel interaction model coupling the g.s. 
of 6 He to the first and second 2 + excitations. This description of the mass-7 systems is 
in many ways alternative to those based on cluster models ( 7 Li as a 3 H+a dicluster, 
etc.) or to those based on more microscopic models, either no-core shell-model or 
Green's function Montecarlo calculations. A comparison amongst different descriptions 
(dicluster model, no-core shell model, and collective-coupling model) for mass-7 nuclei is 
given in Ref. [T7]. In the collective-coupling model, the four A=7 nuclei are described in 
terms of the system nucleon + mass-6-type nucleus with a potential whose parameters 
are listed in Table [U 

7 Li^p + 6 He(0i[g.s.];2i[1.78MeV];2i[5.6MeV]) 

7 He^n + 6 He(0+[g.s.];2+[l.78MeV];2+[5.6MeV]) 

7 Be^n + 6 Be{0+[g.s.];2+[l.70MeV];2+[5.6MeV}) 

7 B^p + 6 Be{0i[g.s.];2i[l.70MeV];2+[5.6MeV}) 



In Table [3] are shown the results obtained .17] with that CC nucleon-nucleus poten- 
tial model. One CC potential model has been used for all four nuclides. The results 



8 



V = 


-36.817 V u = -1.2346 


V U = 


14.9618 Vis = 0.8511 


R = 


2.8 fm; a = 0.88917 fm; 


ft = 


0.7298 



Table 2: Parameter values of the nucleon-nucleus(mass-6) potential. 
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-9.975 


-9.975 


-10.676 


-11.046 












1- 
2 


-9.497 


-9.497 


-10.246 


-10.680 












7 - 
~ 


-5.323 [0.069] 


-5.323 


-6.106 [0.175] 


-6.409 












2" 


-3.371 [0.918] 


-3.371 


-3.946 [1.2] 


-4.497 


J" 


7 Hc 
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5- 
2 


-2.251 [0.08] 


-0.321 


-3.466 [0.4] 


-1.597 




Exp. 


Theory 


Exp. 


Theory 


3- 
2" 


-1.225 [4.712] 


-2.244 






3~ 
2" 


0.445 (0.15) 


0.43 (0.1) 


2.21 (1.4) 


2.10 (0.19) 


1- 
2 


-0.885 [2.752] 


-0.885 




-2.116 


7- 
2 




1.70 (0.03) 




3.01 (0.11) 


7- 
2 


-0.105 [0.437] 


-0.405 


-1.106 [?] 


-1.704 
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1.0 (0.75) ?" 


2.79 (4.1) 




5.40 (7.2) 


3 " 






-0.776 [1.8] 


-3.346 


5- 
2 


3.35 (1.99) 


3.55 (0.2) 




5.35 (0.34) 


3~ 
2 


1.265 (0.26) 


0.704 (0.056) 


0.334 (0.32) 


-0.539 


3- 
2 


6.24 (4.0) ?" 


6.24 (1.9) 






1- 
2 




1.796 (1.57) 




0.727 (0.699) 
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3.7 (0.8) ?" 


2.981 (0.99) 




1.995 (0.231) 
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4.7 (0.7) ?" 


3.046 (0.75) 




2.009 (0.203) 












2 




5.964 (0.23) 




4.904 (0.150) 












7~ 
2 




6.76 (2.24) 


6.5 (6.5) ? a 


5.78 (1.65) 













Table 3: Experimental data and theoretical MCAS results for 7 Li and 7 Be states (left 
table), and for 7 He and 7 B (right table). All energies are in MeV and relate to scattering 
thresholds for nucleon+ 6 He or nucleon+ 6 Be. For states labeled by "? a " spin-parity 
attributes are unknown or uncertain. 

for 7 Li/ 7 Be (and for 7 He/ 7 B) differ solely by the effect of the central Coulomb field. 
Instead, if we compare results for the pairs 7 Li/ 7 He and 7 Be/ 7 B, they differ solely for 
the different action of the OPP term. For the two mass-7 bound systems, Pauli block- 
ing is assumed in the 0si/ 2 shells and all the remaining shells are considered allowed, 
while for the two unbound systems a more complex OPP scenario is considered: the 
0si/2 shells are blocked, the p shells are hindered, and only the higher shells are com- 
pletely allowed. This hindrance of the p shells could reflect the anomalous interaction 
of the neutron with 6 He, which is a neutron halo. A similar situation could occur in 
the mirror case of 7 B, with the interaction of a proton with a 6 Be-type core. In our 
calculation, we have made the hypothesis of a Pauli hindrance in the p-shells defined by 
the following parameters: A(0p 3/2 [0| ; 2+; 2+]) = 17.6MeF, A(0p 1/2 [0^]) = 36.0M e y, 
A(0pi/ 2 [2^; 2^]) = 5.6MeV. The extended nature of the even-even mass-6 subsystems, 
either weakly-bound ( 6 He) or unbound ( 6 Be), is approximately reflected also in the 
geometric size of the potential parameters of Table [2] with rather extended radius, 
diffuseness, and quadrupole deformation. 

4 Conclusions 

The Sturmian-based MCAS approach has been applied to coupled-channel problems at 
low energy, using phenomenological potentials with macroscopic, collective-type, cou- 
plings. But the method is sufficiently flexible that it could be applied also in the 
presence of nonlocal potentials, such as those microscopically generated. Using simpli- 
fied collective-type couplings, we have applied the approach to stable nuclei, as well as 
to weakly-bound and to unbound (with respect to the nucleon emission threshold) light 
nuclei. In the few cases considered, interesting results have been obtained, sometimes 
with rather good reproduction of bound spectra and scattering observables, and often 
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with predictions that could stimulate new experiments. In this approach, the highly 
nonlocal OPP term is crucial in order to include, macroscopically, the effects of Pauli 
exchanges. Finally, for weakly-bound or unbound light systems, the concept of Pauli 
hindrance is suggested. This implies that the nuclei have partially occupied p-shells or 
p-wave proto-shells wherein the Pauli principle neither forbids nor allows occupancy. 
To some extent, it represents a suppression of access to phase-space. 
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